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ABSTRACT 


Several  algorithms  based  on  autocorrelation  matching  of  multiple 
hydrophone  elements  in  a  vertical  line  array  have  been  developed  to  localize  a 
broadband  transient  signal.  An  earlier  developed  frequency-domain 
autocorrelation  matching  (FACM)  algorithm  was  based  on  autocorrelation 
matching  of  only  a  single  hydrophone.  The  success  and  robustness  of  this 
algorithm  in  the  presence  of  environmental  mismatch  was  the  motivation  to  adapt 
it  to  include  the  additional  information  of  multiple  hydrophones.  The  new 
algorithms  developed  were  based  on  joint  autocorrelation  matching,  specifically 
depth-  and  frequency-domain  autocorrelation  matching  (ZFACM),  wavenumber- 
and  frequency-domain  autocorrelation  matching  (KzFACM),  and  an  incoherent 
summation  of  the  FACM  results  of  all  the  elements  in  a  vertical  line  array 
(IFACM).  These  algorithms  were  tested  in  simple,  shallow  water  environments 
with  and  without  mismatch  in  the  specification  of  acoustic  parameters  (e.g., 
bathymetry  and  sound  speed). 

The  results  suggest  that  the  use  of  the  additional  information  from  multiple 
elements  does  improve  both  the  accuracy  and  robustness  of  the  localizations. 

All  of  the  algorithms  produced  similar  results,  although  the  IFACM  and  the 
KzFACM  algorithms  appeared  to  perform  slightly  better  than  the  ZFACM 
algorithm  in  the  presence  of  mismatch.  However,  the  relative  performance  of  the 
algorithms  appeared  to  be  sensitive  to  the  environment  and  placement  of  the 
source  and  receivers  in  the  waveguide. 
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I.  INTRODUCTION 


Today,  matched  field  processing  (MFP)  is  one  of  the  most  commonly  used 
and  successful  research  techniques  for  determining  the  lacking  range  and  depth 
information  for  passive  sonar  systems.  In  MFP,  a  comparison  is  made  between 
some  features  of  a  received  signal  (e.g.,  amplitude  and  phase)  and  those  of 
replica  signals  generated  by  a  numerical  model  predicting  the  field  due  to  a 
source  positioned  at  each  point  in  the  search  grid.  A  match  is  obtained  when  the 
correlation  of  the  two  signals,  synthetic  and  true,  is  at  a  maximum. 

Since  1993,  the  passive  transient  localization  (PTL)  project  at  the  Naval 
Postgraduate  School  (NPS)  has  been  studying  various  algorithms  for  localizing 
transient  sources  based  on  methods  related  to  MFP.  The  interest  in  transient 
signals  is  based  on  the  fact  that  these  signals  can  represent  several  (some  very 
important)  instances  of  a  ship’s  routine,  such  as  opening  of  hatches,  starting  of 
pumps,  popping  of  the  hull,  etc. 

The  primary  acoustic  signal  feature  on  which  most  of  the  NPS  effort  has 
focussed  has  been  the  signal  autocorrelation.  The  motivation  for  this  was  based 
partly  on  the  lack  of  absolute  travel  time  information  as  well  as  the  ability  of  the 
autocorrelation  to  minimize  the  influence  of  noise  at  non-zero  lag  values. 
Originally,  time-domain  autocorrelation  matching  techniques  were  analyzed 
(Miller,  et  al,  1996,  de  Kooter,  1997).  Subsequent  analysis  was  also  performed 
on  frequency-domain  autocorrelation  matching. 

One  of  the  latest  studies  for  the  project  included  the  development  of  MFP 
algorithms  for  a  single  hydrophone  receiver  and  a  broadband,  transient-like  point 
source  (Brune,  1998).  The  performance  of  these  algorithms  was  examined 
numerically  in  environmental  mismatch  situations.  One  of  the  results  of  this 
study  was  that  the  frequency-domain  autocorrelation  matching  (FACM)  algorithm 
was  most  robust  for  environmental  mismatch  and  in  the  presence  of  noise.  This 
was  most  likely  due  to  less  sensitivity  to  phase  mismatch  at  higher  frequencies. 
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Another  recent  study  for  the  PTL  project  examined  the  success  of  the 
FACM  algorithm  on  measured  transient  data  in  a  shallow  water  environment 
(Correa,  1998).  In  this  experiment,  a  multiple  element  vertical  line  array  was 
used.  However,  only  single  element  analysis  was  carried  out.  In  that  thesis,  one 
of  the  recommendations  for  future  work  was  an  adaptation  of  the  FACM 
localization  algorithm  to  include  multiple  element  information.  That  was  the 
motivation  for  the  present  study. 

The  objective  of  this  thesis  is  to  localize  (in  range  and  depth)  a  broadband, 
transient-like  point  source  with  a  vertical  line  array  of  16  hydrophone  receivers. 
Using  the  geometry  of  the  vertical  line  array,  and  re-coding  the  previously 
developed  FACM  algorithm,  a  number  of  new  MFP  algorithms  have  been 
developed.  These  algorithms  are  a  depth-  and  frequency-domain  autocorrelation 
matching  (ZFACM)  algorithm,  a  wavenumber-  and  frequency-domain 
autocorrelation  matching  (KzFACM)  algorithm,  and  an  incoherent  summation  of 
FACM  results  of  the  1 6  receiver  hydrophones  in  the  array.  The  robustness  of  the 
new  algorithms  has  been  tested  numerically  in  different  environments  and  with 
different  environmental  mismatches. 

The  remainder  of  this  thesis  consists  of  four  chapters.  In  Chapter  II,  the 
localization  algorithms  are  described.  In  Chapter  III,  the  propagation  model  that 
was  used  and  the  necessary  data  manipulations  are  discussed.  Chapter  IV 
describes  the  setup  of  the  numerical  experiments  and  the  results.  Finally, 
Chapter  V  provides  conclusions  and  recommendations  for  future  work. 


II.  LOCALIZATION  ALGORITHMS 


As  mentioned  in  the  Introduction,  the  realization  of  the  algorithms  for  this 
thesis  is  an  adaptation  of  previous  work.  Earlier  developed  algorithms,  based  on 
MFP  principles,  dealt  with  a  single  hydrophone  receiver  and  a  transient-like, 
broadband  point  source.  The  matching  took  place  between  the  received  signal 
and  a  replica.  A  replica  is  basically  a  prediction  of  the  field  due  to  a  source 
positioned  at  each  point  in  the  search  grid.  A  numerical  model  that  will  be 
discussed  later  generated  these  predictions.  From  this  previous  work,  the 
frequency-domain  autocorrelation  (FACM)  algorithm  was  concluded  to  give  the 
best  results.  Best  results  were  defined  as  most  robust  for  environmental 
mismatch,  which  was  most  likely  due  to  less  sensitivity  to  phase  mismatch  at 
higher  frequencies.  Therefore,  the  FACM  algorithm  is  best  suited  for  larger 
bandwidths  and  higher  frequencies. 

The  new  algorithms  are  a  continuation  of  this  approach,  taking  the  depth 
information  of  the  multiple  hydrophone  receivers  into  account.  The  transient-like 
point  source  is  assumed  to  be  omnidirectional.  In  order  to  reduce  the  amount  of 
work  to  a  manageable  size,  the  reciprocity  principle  was  used.  This  principle 
states  that  in  an  environment  without  any  time  variations  (e.g.,  currents),  and  with 
identical  densities  at  location  A  and  location  B,  the  acoustic  pressure  at  B  from 
an  omnidirectional  source  at  A  is  identical  to  the  acoustic  pressure  at  A  from  an 
equivalent  source  at  B.  Therefore,  the  replicas  can  be  generated  from 
predictions  of  the  acoustic  propagation  from  the  hydrophone  receiver  location  to 
all  possible  source  locations  in  the  search  grid. 
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A.  FREQUENCY-DOMAIN  AUTOCORRELATION  MATCHING 


Since  autocorrelation  matching  was  initially  developed  for  use  in  the  time- 
domain,  it  is  more  convenient  to  start  the  explanation  in  this  domain  and  convert 
the  results  to  the  frequency-domain.  Suppose  we  have  a  complex  time  series 
P(t)  representing  a  detected  transient  arrival.  The  signal  autocorrelation  in  the 
time-domain  can  be  expressed  as 


TpM=  (2-1) 


or.  in  terms  of  the  frequency-domain  response,  as 


where 


r„(r)=  \p-(f)p(f)e-^*’df . 

— oo 

p{f)=  ]p{ty^’^dt. 


(2-2) 


(2-3) 


The  normalized  autocorrelation  function  is  defined  as 


^pp  (^)  “ 


||p(/f# 


(2-4) 


For  the  predicted  replica  signal  R(t),  a  similar  derivation  leads  to 


(2-5) 


4 


The  time-domain  autocorrelation  matching  (TACM)  algorithm  is  based 
upon  the  inner  product  of  the  two  quantities  in  Eqs.  (2-4)  and  (2-5).  Normalizing 
this  result,  it  can  be  expressed  (in  the  time-domain)  as 


Avica/  ~ 


1/2 


(2-6) 


Note  that  this  technique  attempts  to  match  the  relative  correlations  of 
separate  multipaths  in  the  time-domain  signal.  These  relative  correlations  occur 
at  lag  values  corresponding  to  their  temporal  separation  at  the  receiver.  The 
temporal  separations  changed  with  the  sources  range  and  depth  position  in  the 
environment  and  provided  the  information  necessary  to  uniquely  localize  the 
source.  Previous  analyses  with  algorithms  based  on  this  method  were  found  to 
be  highly  sensitive  to  mismatch  at  high  frequencies.  Specifically,  if  the  relative 
arrival  times  of  different  multipaths  are  off  by  only  a  few  milliseconds,  the  phase 
mismatch  will  be  quite  large  at  frequencies  of  interest  for  transients  (~  1  kHz). 

In  the  frequency-domain,  however,  such  frequency-dependent  effects  do 
not  occur  where  only  relative  frequencies  are  relevant.  Thus,  a  frequency- 
domain  autocorrelation  matching  algorithm  was  developed.  By  analogy,  the 
autocorrelation  matching  in  the  frequency-domain  becomes 


A 


FACM 


X  ' 


(2-7) 
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where  and  C^-  are  the  frequency-domain  autocorrelations  of  the  transient 
signal  and  of  the  replica,  respectively.  Since 


y(tfdi\\p(tfdt 


\dt\R{tf  \dt\P(t'f 


'\R{,fdt]\P{tfdt 


]\R{lfd,]\P(,fdt 


(2-8) 


Eq.  (2-7)  can  be  rewritten  as 


J|i?(r)|^|p(r|' 


dt 


^FACM 


X 


(2-9) 


To  reduce  the  influence  of  noise  and  to  emphasize  the  comparison 
between  the  cross  terms,  it  is  necessary  to  remove  the  zero-lag  of  the 
autocorrelation  function.  This  is  equivalent  to  removing  the  mean  squared 
amplitude  of  the  signal  in  the  time-domain.  Therefore,  defining 


where  in  both  cases  the  average  is  over  all  times,  Eq.  (2-9)  becomes 


A 


FACM  ~ 


|>2 


(2-12) 


In  order  to  find  the  optimal  match,  a  search  parameter  t  will  be 
introduced  which  allows  us  to  slide  the  replica  signal  in  time.  This  procedure  is 
necessary  since  absolute  time  information  is  not  available.  The  FACM  algorithm 
is  then  based  on  evaluation  of  the  function 


^FACM  ~  T 


f|R'('+^ 

y{>u 

[.(I 

R'itXd,] 

(2-13) 


where  max^  is  the  maximum  value  of  the  function  for  any  value  of  the  search 
parameter. 


In  previous  thesis  work,  a  MATLAB®  implementation  of  the  algorithm 
based  on  Eq.  (2-13)  was  used  to  generate  the  ambiguity  surfaces  for  a  single 
hydrophone  receiver.  The  sliding-  t  operation  was  performed  by  a  multiplication 
in  the  frequency-domain  and  an  inverse  FFT  operation,  as  will  be  shown  by  the 

following  derivation.  Let  g(r)  =  and  h{t)  =  |P'(r)|^ ,  then 


and 


g{l)=  \e-'^*G(f)df  , 

—CO 

eo  eo 

h{t  +  r)  =  J H{f)df  =  J H{f)e-‘^’^^df  . 


The  numerator  of  Eq.  (2-13)  becomes 


oo 

r(r)=  ^h(t  +  T)g(t)dt 


Since  h{t)  is  a  real  function,  then  h{t)-h'{t)  and 


i{t)=  jh"{t  +  T)g(t)dt 


\dt 


df'df 


(2-14) 


(2-15) 


(2-16) 
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(2-17) 


—oo 

Thus,  after  computing  corresponding  frequency- 

domain  representations  g(/)  and  //(/),  the  function  l(r)  can  easily  be 
computed.  The  maximum  value  of  /(r)  then  corresponds  to  the  maximum  of  the 
numerator  in  Eq.  (2-13). 

B.  INCOHERENT  SUMMATION 

The  vertical  line  array  has  16  omni-directional  hydrophone  receivers.  The 
most  straightforward  extension  of  the  FACM  algorithm  is  the  incoherent 
summation  of  the  FACM  results  of  each  of  the  single  hydrophone  receivers  of  the 
array.  As  mentioned  earlier,  the  FACM  results  are  displayed  as  ambiguity 
surfaces.  Therefore,  the  incoherent  summation  (IFACM)  algorithm  can  best  be 
seen  as  the  averaged  overlay  result  of  16  ambiguity  surfaces  of  single 
hydrophone  receivers.  This  can  be  defined  simply  as 


— oo 
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C.  DEPTH-  AND  FREQUENCY-DOMAIN  AUTOCORRELATION 


MATCHING 

A  more  advanced  extension  of  the  FACM  algorithm  is  the  depth-  and 
frequency-domain  autocorrelation  matching  (ZFACM)  algorithm,  where  the 
hydrophone  receiver  depths  are  incorporated  into  the  autocorrelation  calculations 
as  a  second  ‘dimension’.  This  can  also  be  thought  of  as  joint  autocorrelation 
matching. 

The  following  derivation  of  the  ZFACM  algorithm  has  the  same  flow  as  the 
‘one-dimensional’  FACM.  This  time,  the  derivation  will  start  directly  in  the  depth- 
and  frequency-domain.  A  necessary  additional  remark  beforehand  concerns 
Fourier-pairs.  The  depth-  and  vertical  wavenumber-values,  ,  form  a  Fourier- 
pair,  similar  to  the  probably  better  known  Fourier-pair,  frequency-  and  time- 
values.  Thus,  we  will  simplify  our  development  by  taking  advantage  of  such 
symmetries. 

Because  all  the  hydrophone  receivers  are  located  at  a  specific  depth,  the 
transient  arrival,  as  detected  by  the  whole  vertical  line  array,  can  be  expressed 
by  a  complex  time-  and  depth-series  p{t,z).  With  a  Fourier  transform  on  the 
time-domain  component  of  this  series,  it  will  become  a  complex  depth-  and 
frequency-series  P{z,f). 

The  signal  autocorrelation  in  the  depth-  and  frequency-domain  may  be 

written 


Fpp{^,<l>)=  J  \P*{z,f)P{z^^,f  +  (l>)lzdf  .  (2-19) 
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Defining 


j  oa  oo 

P{z,  /)  =  ^  J  J  (2-20) 

—00—00 

P{z,f)  =  F-‘[p{k^,t)] 
and 

P(k,,t)^  j  (2-21) 

—00—00 

/>(*„<)  =  F’[/>(z,/)]  , 

where  F^and  F'^  indicate  the  two-dimensional  forward  and  inverse  Fourier 
transforms,  respectively,  Eq.  (2-19)  can  be  rewitten  as 


dzdf 


1 

2n 


J  j  J  ^dtdt'dk.dk[P^(k.j)p[k[,t'y^”'^'e'^'^ 


1 

271 


=  -^  J  J P'(k.^,t)p[k.^,ty^’^^e‘^-~dtdk.^ 


(2-22) 
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The  normalized  joint  autocorrelation  becomes 


Cpp{^,0)  = 


±]]\p{k,.,fdk,d, 


(2-23) 


For  the  predicted  replica  signal  R{z,f),  a  similar  derivation  leads  to 


C..  (#,«*) 


(2-24) 


Since  we  want  to  emphasize  the  matching  of  the  non-zero  lag  values, 
reducing  the  influence  of  noise,  these  zero-lag  values  should  be  removed.  This 
can  be  accomplished  by  removing  the  means  from  the  wavenumber-  and  time- 
domain,  i.e. 


and 


J_ 

In 


J  dk.dt 


In 


J  dk,dt 


(2-25) 


(2-26) 
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where 


and 


(2-27) 


(2-28) 


In  both  cases  the  averages  are  taken  overtime  and  wavenumber.  The  ambiguity 
surface  is  then  based  on 


'■ZFACM 


K  ■ 


(2-29) 


where 


oo  oo 


j  i  j  J\R'{k..tf<!-‘^'‘e-‘‘’*dk^dt  j  J\p'{kUi'f  e‘^'‘e‘^'dk[dt' 


In 


tA  j  ]\K'(KA^dk.df]  ]\p'(k  .Ipk^d, 

y^/i-  J  _oo— oo 


-^1  j  j|p'(t  .tjdkji 

\^n )  _„_oo 
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{iTlf 


CO  oo 


1  ‘ 

iTt  • 

o  oo 

fj|4 

00—00 

■A 

dk.dt 

j  ]\R’{K.‘ 

]\A  A 

—00—00  —00—00 


dk,dt 


Then,  also 


—00—00 


and 


1 

2n: 


(2-30) 


(2-31) 


(2-32) 


Finally,  the  ZFACM  function  can  be  expressed  as 

^ZFACM  ~ 


y(k„t 

1  dkMt 

Xdk,dt 

'y}\AA 

2 

dk^dt 

^  J  ]\Ak-Adk,dt  j-  J  ]\p'(k,jXdkJt 


}  —00—00  —00—00 


14 


iX 


(2-33) 


dk.dt 


J  ^\R'(k,  ,tXdkJt- ^  ^\p\k.  ,tX  dkjt 


Like  with  the  FACM  derivation,  since  absolute  time  information  is  not 
available,  a  search  parameter  r  will  be  introduced  in  order  to  find  the  optimal 
match  in  time.  Also,  a  search  parameter  K  will  be  introduced  as  the 
wavenumber-domain  equivalent.  Since  the  absolute  depths  of  the  receivers 
(hence  absolute  wavenumber)  are  exactly  known,  the  second  search  parameter 
seems  to  be  obsolete.  However,  the  autocorrelation  matching  has  become  a 
‘two-dimensional’  matching,  therefore  both  search  parameters  are  introduced. 
These  search  parameters  allow  us  to  slide  the  replica  signal  in  wavenumber  and 
time.  The  ZFACM  function  then  becomes 


^ZFACM 


J  '!*<*  =  <*  j  f\P'(kY' 


1X2 


(2-34) 


The  sliding  operation  is  again  performed  in  the  alternate  domain  in  a  manner 
similar  to  that  described  in  the  previous  description  of  the  FACM  algorithm. 
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D.  WAVENUMBER-  AND  FREQUENCY-DOMAIN 


AUTOCORRELATION  MATCHING 

Another  advanced  extension  of  the  FACM  algorithm  is  the  wavenumber- 
and  frequency-domain  autocorrelation  matching  (KzFACM)  algorithm.  This  time, 
the  vertical  wavenumber-values  are  incorporated  into  the  autocorrelation 
calculations  as  the  second  ‘dimension’. 

The  derivation  of  the  KzFACM  algorithm  has  exactly  the  same  flow  as  the 
ZFACM  derivation,  and  therefore  the  larger  part  of  it  will  be  left  out.  A  complex 
time-  and  depth-series  P{t,  z) ,  as  detected  by  the  whole  vertical  line  array,  can 
express  the  transient  arrival.  With  a  Fourier  transform  on  both  domain 
components  of  this  series,  it  will  become  a  complex  wavenumber-  and 
frequency-series  p{k.,f). 

A  similar  derivation  as  for  the  ZFACM  function  leads  to  the  final  result  for 
the  KzFACM  function. 


^KzFACM  ~ 


R'{z  +  ^,  r  +  t)|  ^  |f'(z,  r)|  ^  dzdt 


J  ^\R'{z,t\ dzdt  ■  J  ^\p'{z,t\ dzdt 


(2-35) 


where  ^  and  r  are  the  search  parameters  for  the  depth-  and  time-domain, 
respectively.  Again,  since  the  autocorrelation  matching  has  become  a  ‘two- 
dimensional’  matching,  both  search  parameters  are  introduced.  These  search 
parameters  allow  us  to  slide  the  replica  signal  in  depth  and  time. 
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It  is  important  to  note  the  physical  interpretation  of  these  autocorrelation 
matching  algorithms.  As  has  been  shown,  the  autocorrelation  matching  in  one 
domain  may  be  computed  from  the  square  amplitude  in  the  Fourier  domain. 
Thus,  for  example,  the  FACM  algorithm  is  actually  performing  a  squared 
amplitude  (envelope)  matching  in  the  time-domain.  In  the  context  of  such 
envelope  matching,  it  is  perhaps  more  obvious  that  this  algorithm  will  be  less 
sensitive  to  time-domain  phase  mismatch.  Therefore,  an  examination  of  the 
previous  work  suggests  that  the  time-domain  envelopes  provide  a  more  stable 
acoustic  signature  for  transient  localization  than  the  corresponding  frequency- 
domain  envelopes. 

The  relationship  between  domains  should  be  kept  in  mind  when 
evaluating  the  following  results  presented  in  this  thesis.  Specifically,  the  depth- 
and  frequency-domain  autocorrelation  matching  is  equivalent  to  a  squared 
amplitude  matching  in  the  wavenumber-  and  time-domain.  Furthermore,  it  is 
worthwhile  noting  that  the  results  in  the  wavenumber-domain  may  also  be 
associated  with  similar  algorithms  employed  in  beam-space. 
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III.  NUMERICAL  METHOD 


A.  MONTEREY-MIAMI  PARABOLIC  EQUATION  (MMPE) 
PROPAGATION  MODEL 

In  this  chapter,  the  propagation  model  used  to  predict  the  arrival  structure 
replicas  at  a  receiver  location  due  to  a  point  source  at  a  test  location  will  be 
described.  Also,  the  methods  for  processing  the  outputs  of  the  model  to  produce 
the  time-domain  arrival  structure  predictions  will  be  explained. 

In  order  to  generate  the  replicas  that  simulated  the  received  signals  from  a 
synthetic  source  on  defined  points  of  a  two  dimensional  search  space  (depth 
versus  range),  the  Monterey-Miami  Parabolic  Equation  (MMPE)  propagation 
model  (Smith,  1999)  was  used.  In  this  chapter,  the  general  theory  behind  the 
parabolic  equation  model  will  be  introduced,  as  well  as  the  method  used  for  its 
implementation  -  the  split-step  Fourier  method  (Hardin  and  Tappert,  1973). 

The  inhomogeneous  wave  equation  for  the  acoustic  pressure  p{x,t)  in  a 
medium  with  sound  speed  c(^)  and  density /o{jc)  can  be  expressed  as  (Jensen, 
etal,  1994) 


1  d^p{x,t) 
c(3c)'  dt' 


(3-1) 


The  parabolic  equation  model  is  based  on  an  approximation  of  the  Helmholtz 
wave  equation  in  a  cylindrical  coordinate  system.  Because  of  the  ocean’s 
relative  shallowness  compared  to  horizontal  propagation  distance  for  the  majority 
of  the  environments,  it  is  well  suited  for  a  description  in  cylindrical  coordinates. 
Assuming  a  time  harmonic  solution,  the  Helmholtz  equation  takes  the  form 
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1  d  f^  dp(r,z,^,f))  1  d^p[r,z,<l>,f)  ^d'^p[r,z,(pj) 


dr 


dr  J  r- 


^  d^^  dz^ 

+  kln^ (r, z, (l>)p(r, z,^,f)=  -4;tPo j(? - rj  , 


(3-2) 


where 

P{r,z,^,t)=  p{r,z,(p,f)e~‘^'^  .  (3-3) 

The  reference  wavenumber  is  related  to  a  reference  sound  speed,  ,  by 


and  the  acoustic  index  of  refraction  is  defined  by 


(3-4) 


n  = 


(3-5) 


Note  that  in  this  derivation  the  density  variations  are  neglected.  By  defining  the 
effective  index  of  refraction  which  contains  the  appropriate  additional  terms,  the 
influence  of  the  density  differences  at  the  water-bottom  interface  can  be  included 
(e.g.,  Tappert,  1977). 

By  assuming  that  the  ocean  acts  as  a  waveguide  with  a  cylindrical 
coordinate  system,  acoustic  energy  is  mainly  propagated  outward  from  a  source 
in  the  horizontal  direction.  Therefore,  the  pressure  field  can  be  approximated  by 

p{r,z,(p,f)  =  ^(r,z,0,f)H^'^{kor)  ,  (3-6) 
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where  is  a  slowly  varying  envelope  modulating  the  outgoing  zero-th 

order  Hankel  function  of  the  first  kind  Taking  advantage  of  the  far-field 

{k^r  » l)  asymptotic  approximation  of  the  Hankel  function,  Eq.  (3-6)  can  be 
rewritten  as 


(3-7) 


normalized  such  that  at  the  reference  range  r  =  R^,  |p|  =  Pq  ■  Substituting  Eq. 
(3-7)  into  Eq.  (3-2)  and  dropping  the  source  term  on  the  right  hand  side  gives 


a'p  1  a^T 

8r=  *  ‘ 


a"p 


1 

4^2 


^  =  0  . 


(3-8) 


Neglecting  the  azimuthal  coupling  and  the  far-field  terms,  and  dropping  the  first 
term  due  to  the  slow  modulations  of  the  envelope  function,  Eq.  (3-8)  can  be 
rewritten  as 


/  ay 

kg  dr 


1  a-^ 

2k,  dz^ 


2 


Defining  the  operators 


T  =- 
""  2k 


1  ( a 


2  ^ 


dz 


\  J 


and 


(3-9) 


(3-10) 


(3-11) 
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Eq.  (3-9)  can  be  written  as 


~  =  .  (3-12) 

K  or 

Equation  (3-12)  is  known  as  the  “standard”  parabolic  equation  (SPE)  (Tappert, 
1977),  with  accurate  solutions  limited  to  a  half  beam  width  of  10°  to  20°  for  the 
propagation  angle.  For  this  work,  in  order  to  extend  this  limit  to  40°  to  70°,  a 
higher  order  wide-angle  parabolic  equation  (WARE)  approximation  (Thomson 
and  Chapman,  1983)  is  used.  Its  operators  are  defined  by 


and 


'^WAPE 


=  1- 


(3-13) 


^WAPE  ~  l)  • 


(3-14) 


The  MMPE  uses  the  split-step  Fourier  (SSF)  method  in  order  to  solve  the 
parabolic  equation  numerically.  This  algorithm  integrates  the  solution  in  range  by 
applying  the  and  the  operators  in  the  A:,  -  domain  and  the  z  - 
domain,  respectively,  where  each  operator  is  simply  a  scalar  multiplier.  In  the 
(t,  -  space  the  wide  angle  f^^p^  operator  is  defined  as 


^WAPe{^z)  ^ 


(3-15) 
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The  algorithm  for  stepping  in  range  from  r  to  r  +  Ar  is  the  o[Ar^)  accurate 
centered-step  scheme  (Jensen,  et  al.,  1994),  which  can  now  be  expressed  as 


^(r  +  Ar,  z)  =  -1 


-iX'oAr fwylW;  )p 


(3-16) 


B.  DATA  MANIPULATIONS 

The  output  of  the  MM  PE  model  is  in  the  form  of  the  field  function  'P  (both 
magnitude  and  phase)  and  has  been  referenced  to  unit  magnitude  at  1  m.  The 
field  function  is  computed  at  the  defined  spatial  grid  points.  The  pressure  is 
defined  in  terms  of  the  field  function  'F  by  Eq.  (3-7). 

Broadband  results  were  obtained  by  running  the  MMPE  model  for  all 
discrete  frequencies  in  the  chosen  bandwidth.  The  travel  time  results  were 
realized  by  performing  a  Fourier  synthesis  and  a  necessary  multiplication  by 
some  window-function,  s[f).  The  complex  arrival  structure  of  the  pressure  field 
can  then  be  written  as 


oo 

P[r,z,t)=  ^s{f)p[r,zjY'^’^df  = 

—oo 


^0^0 

4~r 


(3-17) 


By  defining  the  reduced  time  T  =  t- 
neglected,  such  that 


^  r  ^ 


i27f- 


,  the  phase  factor  =e  can  be 


P(r,zj)  =  ^]s(f)v(r,z,(l>y^’^df  . 

yf  -oo 


(3-18) 
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Note  that  the  use  of  a  reduced  time  does  not  influence  the  autocorrelation 
function. 

For  all  algorithms,  the  same  window  function  was  used.  The  modeled 
pulse  transmission  had  a  center  frequency  of  1400  Hz  with  a  bandwidth  of  900 
Hz.  The  temporal  spread  at  maximum  range  was  ~  0.25  seconds.  In  order  to 
avoid  any  errors  due  to  aliasing,  the  minimum  width  of  the  time-domain  window 
should  be  larger  than  this.  Therefore,  512  frequency  bins  were  computed  over 
the  bandwidth,  which  produced  a  time  window  of  0.569  seconds.  Over  this 
bandwidth,  a  simple  Hanning  window  was  applied. 

To  remove  the  mean  from  the  signals,  as  required  by  Eqs.  (2.10)  and 
(2.11),  without  introducing  a  discontinuous  removal  of  the  d.c.-component  in  the 
Fourier  domain,  the  following  approach  was  taken.  To  remove  the  mean  from 

|p(r]j^ ,  for  example,  the  series  |p(r]|^  was  first  transferred  to  the  frequency  - 

domain.  A  smooth  Hanning-taper  was  then  applied  to  the  16  bins  around  d.c., 
followed  by  a  transformation  back  to  the  time-domain.  The  consequence  of  this 
filtering  process  was  to  remove  the  mean  and  filter  only  the  lowest  frequency 
components. 

For  both  the  source  data  and  replica  data,  the  field  information  from  the 
output  of  the  MMPE  propagation  model  was  used  as  input  for  the  autocorrelation 
matching  algorithms.  The  structure  of  the  MMPE  output  file  is  as  follows.  After  a 
header  with  general  information,  the  amplitude  and  phase  values  of  the  field 
function  are  stored  at  every  depth  point  for  the  first  range  step  and  frequency  bin. 
This  is  followed  by  the  same  data  for  the  next  range  step,  and  continued  for  all 
range  steps  out  to  the  maximum  range  for  the  first  frequency  bin.  After  reaching 
the  maximum  range,  the  procedure  begins  again  for  the  next  frequency  bin.  This 
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process  repeats  itself  until  all  frequencies  have  been  computed.  The  result  of 
this  ordering,  from  innermost  to  outermost  loops,  is  depth  -  range  -  frequency. 

In  order  to  develop  the  autocorrelation  matching  algorithms  as  time- 
efficient  as  possible,  the  input  files  should  be  read  in  such  a  way  to  obtain  the  full 
frequency  bandwidth  at  each  spatial  gridpoint.  Therefore,  the  MM  PE  output  files 
had  to  be  reordered  first.  For  this  reordering  procedure  an  existing  FORTRAN® 
algorithm  was  used.  The  result  of  this  reordering  process  produced  data  files 
organized  from  innermost  to  outermost  as  frequency  -  depth  -  range. 
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IV.  EXPERIMENTAL  SETUP  AND  RESULTS 


Several  numerical  experiments  were  conducted  in  order  to  gain  insight 
into  the  performance  of  the  autocorrelation  matching  algorithms.  For  all 
experiments  both  the  transient  and  the  replica  field  functions  were  calculated  with 
the  MMPE  propagation  model.  The  field  function  of  a  replica  represents  those 
values  received  from  a  test  source  positioned  at  each  defined  point  in  the  search 
area.  The  results  of  the  autocorrelation  matching  algorithms  are  presented  in  an 
ambiguity  surface.  For  the  theoretical  situation  where  the  field  functions  of  the 
transient  and  the  replica  are  100%  the  same,  the  ambiguity  surface  will  show  a 
maximum  exactly  at  the  source  location  of  unit  amplitude  (with  proper 
normalization). 

Three  environments  were  selected  to  test  the  robustness  of  the 
algorithms.  The  first  environment  was  mainly  used  to  baseline  the  performance 
of  the  algorithms.  For  this  case,  the  environment  was  exactly  known  for  both  the 
transient  and  the  replicas.  For  each  of  the  other  two  environments,  three 
different  situations  were  tested:  the  all-known  situation,  in  which  the 
environments  for  both  the  transient  and  replica  data  were  exactly  the  same;  an 
environmental  mismatch  in  the  sound  speed  profile;  and  an  environmental 
mismatch  in  the  bottom  bathymetry.  Mismatch  details  will  be  discussed  later. 

The  search  area  had  a  total  range  of  10  km  and  a  depth  of  100  m.  The 
resolution  of  the  ambiguity  surface  is  40  m  in  range  and  1 .56  m  in  depth.  In 
some  situations  a  ‘zoomed-in’  ambiguity  surface  was  calculated  with  a  window  of 
500  m  in  range  and  100  m  in  depth  around  the  source  location.  This  ‘zoomed-in’ 
version  had  a  resolution  of  3.33  m  in  range  and  0.78  m  in  depth. 

For  ail  experiments  the  transient  had  a  1400  Hz  center  frequency  with  a 
900  Hz  bandwidth.  This  bandwidth  was  sampled  over  512  points,  corresponding 
to  a  frequency  resolution  of  1 .7578125  Hz  and  a  time  window  of  width  0.569 
seconds. 
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The  modeled  vertical  line  array  had  a  total  length  of  2.93  m.  This  is 
approximately  the  length  of  a  modern  operational  array.  The  array  consisted  of 
16  omni-directional  hydrophone  receivers.  This  number  was  chosen  in  part 
because  the  calculation  time  of  the  required  Fourier  transforms  can  be 
significantly  reduced  if  the  total  number  of  elements  (in  this  case,  hydrophones) 
is  a  power  of  two.  The  element  spacing  was  0.195  m.  This  element  spacing 
coincided  with  the  nominal  model  depth  grid  spacing,  so  that  all  hydrophones  in 
the  array  could  be  positioned  exactly  on  a  grid  point.  For  each  environment  the 
line  array  was  positioned  at  two  different  depth-ranges,  realizing  a  shallow 
aperture  spanning  19.43  m  to  22.36  m,  and  a  deep  aperture  spanning  54.79  m  to 
57.71  m. 

As  stated  earlier,  the  ambiguity  surfaces  presented  are  the  numerical 
result  of  an  autocorrelation  matching  function,  its  maximum  value  is  unity  since 
the  autocorrelation  matching  function  is  normalized.  This  maximum  value  will 
occur  at  the  actual  source  location,  but  only  in  the  theoretical  situation  that  the 
predicted  field  function  is  identical  to  the  source  field  function.  For  all  the  other 
situations,  more  areas  with  high  correlation  values  can  be  expected  on  the 
ambiguity  surface.  These  areas  are  commonly  referred  to  as  ‘sidelobes’  and 
would  correspond  to  false  targets  in  an  operational  system.  Also,  predicted 
source  positions  are  not  pinpoint  locations  but  finite  area  peaks  of  high 
correlation.  These  areas  give  size  and  shape  to  the  localization  ‘footprint’ 
(Tolstoy,  1993).  The  ambiguity  surface  of  an  algorithm  is  described  by  the 
number  and  magnitude  of  the  sidelobes  and  the  area  of  the  footprint. 

For  the  ambiguity  surfaces,  the  axes  are  defined  as  followed.  Along  the  x- 
axis  the  range  in  km  is  displayed.  This  range  can  either  be  the  full  range  search 
area  of  10  km,  or  the  ‘zoomed-in’  search  area  of  500  m  around  the  source 
location.  Along  the  y-axis  the  depth  in  m  of  the  search  area  is  displayed.  In  all 
cases,  the  depth  ranges  from  0  m  at  the  surface  to  the  maximum  depth  of  100  m. 
The  vertical  line  array  is  not  shown  in  the  ambiguity  surface.  The  array  is  always 
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located  at  a  range  of  0  km,  and  its  depth  is  either  the  shallow  location  or  the  deep 
location.  Next  to  the  ambiguity  surface  a  colorbar  is  shown.  This  colorbar  has  a 
high  limit  equal  to  the  maximum  value  of  the  autocorrelation  matching  function, 
and  a  defined  low  limit.  In  most  cases  this  lower  limit  is  set  to  0.6,  although  in 
some  cases  this  limit  is  adjusted  to  improve  the  readability  of  the  ambiguity 
surface. 


A.  CASE1:  ISOSPEED  SOUND  SPEED  PROFILE 

For  this  case  the  source  was  located  at  a  range  of  5.2  km  from  the  vertical 
line  array  and  at  a  depth  of  33  m.  The  sound  speed  was  range-independent  with 
an  isospeed  profile  of  1 500  m/s.  The  density  of  water  was  assumed  to  be  1 .0 
g/cm^.  The  bottom  was  chosen  to  be  completely  flat  at  a  depth  of  100  m.  Its 
compressional  sound  speed  was  1600  m/s,  the  density  1.5  g/cm^,  and  the 
attenuation  0.1  dB/km/Hz.  This  simple  type  of  environment  is  commonly  referred 
to  as  a  Pekeris  waveguide. 

In  the  following  series  of  figures,  the  ambiguity  surfaces  are  shown  for  the 
shallow  array  for  the  three  different  algorithms:  the  incoherent  summation 
(IFACM)  algorithm  (Fig.  4.1),  the  wavenumber-  and  frequency-domain 
autocorrelation  matching  (KzFACM)  algorithm  (Fig.  4.2),  and  the  depth-  and 
frequency-domain  autocorrelation  matching  (ZFACM)  algorithm  (Fig.  4.3).  For  all 
three,  both  the  full  range  and  the  ‘zoomed-in’  surfaces  are  displayed. 

For  the  IFACM,  the  full  range  ambiguity  surface  is  shown  over  its  full 
dynamic  range  (Fig.  4.1(a))  and  over  an  adjusted  dynamic  range  (Fig.  4.1(b)). 
Since  adjusting  the  dynamic  range  (without  changing  any  of  the  autocorrelation 
matching  results)  will  improve  the  readability,  all  other  ambiguity  surfaces  are 
only  presented  with  adjusted  dynamic  ranges. 
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Figure  4.1 :  For  Case  1 ,  the  ambiguity  surfaces  for  the  IFACM  algorithm  for  the  shallow  array  in  an 
all-known  environment:  (a)  the  full  range  search  area  and  the  full  dynamic  range;  (b)  the  full 
range  search  area  with  an  adjusted  dynamic  range;  (c)  the  ‘zoomed-in’  ambiguity  surface  with  an 
adjusted  dynamic  range. 
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Figure  4.2:  For  Case  1 ,  the  ambiguity  surfaces  for  the  KzFACM  algorithm  for  the  shallow  array  in 
an  all-known  environment;  (a)  the  full  range  search  area;  (b)  the  ‘zoome(j-in’  ambiguity  surface. 
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Figure  4.3:  For  Case  1,  the  ambiguity  surfaces  for  the  ZFACM  algorithm  for  the  shallow  array  in 
an  all-known  environment:  (a)  the  full  range  search  area;  (b)  the  ‘zoomed-in’  ambiguity  surface. 
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Observing  Figs.  4.1  through  4.3  the  following  remarks  can  be  made.  For 
all  three  algorithms  the  maximum  value  of  the  autocorrelation  matching  function 
occurs  at  the  correct  source  location:  5.2  km  and  33  m  deep.  Comparing  the  full 
range  search  areas,  the  IFACM  and  KzFACM  results  show  a  very  good  match 
and  also  a  close  resemblance.  The  ZFACM  result  shows  a  good  match  with 
stronger  sidelobes.  The  footprints  for  IFACM  and  KzFACM  are  almost  equal  in 
shape  and  size,  where  ZFACM’s  footprint  is  larger  in  depth  and  range. 

Comparing  the  ‘zoomed-in’  versions  shows  a  similar  outcome.  The 
IFACM  and  KzFACM  results  have  a  close  resemblance  and  on  both  ambiguity 
surfaces  almost  no  sidelobes  are  present.  Again,  their  footprints  are  practically 
the  same.  The  ZFACM  result  has  more  sidelobes  of  medium-high  magnitude  in 
the  vicinity  of  the  actual  source  location.  Also,  ZFACM’s  footprint  is  smaller  in 
range  but  larger  in  depth. 

In  Fig.  4.4  the  ambiguity  surfaces  are  shown  for  the  deep  array.  For  all 
three  algorithms  the  ‘zoomed-in’  surfaces  with  adjusted  dynamic  range  are 
displayed.  With  regards  to  this  figure  the  following  remarks  can  be  made.  For  all 
three  algorithms  the  maximum  value  of  the  autocorrelation  matching  function  is 
at  the  correct  source  location.  Again,  the  IFACM  and  KzFACM  results  show  a 
very  good  match  and  also  a  close  resemblance.  No  other  sidelobes  are  present 
on  these  ambiguity  surfaces.  The  ZFACM  result  is  a  good  match  with  only  one 
sidelobe  of  medium-high  magnitude  directly  under  the  source  location.  The 
footprints  for  IFACM  and  KzFACM  are  nearly  the  same,  where  ZFACM’s  footprint 
is  marginally  larger  in  depth. 

Comparing  these  ‘zoomed-in’  versions  for  the  deep  array  with  the 
‘zoomed-in’  ambiguity  surfaces  for  the  shallow  array  (Figs.  4.1(c),  4.2(b)  and 
4.3(b)),  the  following  can  be  observed.  For  all  three  algorithms,  the  match  for  the 
deep  array  is  better  than  for  the  shallow  array,  since  fewer  sidelobes  of  medium 
magnitudes  are  present.  Especially  the  ZFACM  match  has  improved.  Also  for 
all  algorithms  the  deep  array  footprints  are  (marginally)  smaller  in  range  and 
twice  as  small  in  depth. 
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Figure  4.4:  For  Case  1 ,  the  ‘zoomed-in’  ambiguity  surfaces 
environment:  (a)  the  IFACM  algorithm;  (b)  the  KzFACM  algoril 
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This  performance  difference  between  the  two  arrays  may  be  explained  as 
follows.  For  an  isospeed  environment,  an  array  closest  to  the  middle  of  the  water 
column  will  intercept  more  multipath  structure  than  an  array  further  away  from  the 
middle.  In  this  experiment,  the  water  column  is  100  m  deep.  The  deep  array  is 
roughly  located  at  a  depth  of  55  m,  which  is  nearly  halfway  in  the  water  column. 

The  shallow  array,  on  the  other  hand,  is  roughly  located  at  a  depth  of  20 
m,  which  is  approximately  at  1/5  of  the  water  column.  The  central  position  in  the 
water  column  of  the  deep  array  will  make  more  interception  of  multipath  structure 
possible,  resulting  in  a  better  performance  than  the  shallow  array. 

Recalling  the  fact  that  Case  1  is  the  theoretical  case  where  transient  and 
replica  signals  are  100%  the  same,  one  would  expect  the  maximum  value  of  the 
autocorrelation  matching  algorithms  to  be  unity.  However,  this  is  not  the  case,  as 
can  be  seen  in  these  first  ambiguity  surfaces.  The  most  likely  factor  causing  this 
is  the  discrepancy  between  the  source  location  and  the  grid  point  locations.  The 
setup  of  the  experiment  is  chosen  in  such  a  way  that  all  receiver  locations,  range 
and  depth,  exactly  coincide  with  the  grid  point  locations.  For  the  source 
locations,  the  source  range  usually  coincides  with  grid  point  range  locations. 
However,  the  source  depth  seldom  coincides  with  the  grid  point  depth  locations. 
This  small  source  location  error  will  contribute  to  the  autocorrelation  value  error. 
Also,  at  a  smaller  scale,  round-off  errors  will  contribute  to  the  autocorrelation 
value  error. 


B.  CASE  2:  POSITIVE  SSP  GRADIENT 

In  Case  2,  the  source  was  located  at  a  range  of  4.6  km  from  the  vertical 
line  array  and  at  a  depth  of  27  m.  The  sound  speed  was  range-independent  with 
a  positive,  linear  gradient  (upward  refracting)  sound  speed  profile  of  1497  m/s  at 
the  surface  and  1499  m/s  at  a  depth  of  100  m.  The  density  of  water  was 
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assumed  to  be  1 .0  g/cm^.  The  bottom  was  chosen  to  be  a  somewhat  ‘bumpy’ 
bottom  with  the  following  depths:  100  m  at  0  km;  102  m  at  2.5  km;  99  m  at  5  km; 
100  m  at  7.5  km;  and  102  m  at  10  km.  Its  compressional  sound  speed  was  1595 
m/s,  the  density  1.52  g/cm^,  and  the  attenuation  0.1 1  dB/km/Hz. 

In  this  case,  three  different  environmental  situations  are  considered.  For 
the  first  situation,  the  environment  is  completely  known  and  the  transient  and 
replica  field  functions,  determined  by  the  MMPE  propagation  model,  have  exactly 
the  same  input  parameters.  In  other  words,  there  is  no  environmental  mismatch. 
Figure  4.5  shows  the  ambiguity  surfaces  for  the  shallow  array  for  all  three 
algorithms. 

Observing  this  figure  the  following  remarks  can  be  made.  For  all  three 
algorithms  the  maximum  value  of  the  autocorrelation  matching  function  occurs  at 
the  correct  source  location;  4.6  km  and  27  m  deep.  All  three  algorithms  show  a 
close  resemblance,  since  the  sidelobes  occur  at  exactly  the  same  locations.  The 
difference  between  the  results  is  mostly  evident  in  the  magnitude  of  the  sidelobes 
and  marginally  in  the  peak  values.  The  IFACM  and  KzFACM  results  show  a 
good  match,  since  only  sidelobes  of  medium-low  magnitudes  occur  on  the 
ambiguity  surface.  The  ZFACM  result  is  a  moderate  match,  because  more 
sidelobes  of  medium-high  magnitudes  are  displayed.  For  all  three  algorithms 
these  sidelobes  mainly  occur  behind  the  source  location  (seen  from  the  vertical 
line  array),  acting  like  ‘shadows’.  The  footprints  are  the  same  in  range  for  all 
three  algorithms,  although  ZFACM’s  footprint  is  marginally  larger  in  depth. 
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Figure  4.5:  For  Case  2,  the  ambiguity  surfaces  for  the  shallow  array  with  no  environmental 
mismatch:  (a)  the  IFACM  algorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM  algorithm  results. 
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Comparing  the  shallow  array  results  of  Case  2  with  the  shallow  array 
results  of  Case  1  (Figs.  4.1(b),  4.2(a)  and  4.3(a)),  the  following  can  be  observed. 
Although  in  each  case  the  transient  and  the  replica  signals  were  100%  the  same, 
an  isospeed  environment  (Case  1 )  seems  to  provide  a  better  match  for  all  three 
algorithms  than  a  positive  SSP  environment  (Case  2).  The  footprint  sizes  are 
approximately  the  same.  However,  in  Case  2  more  sidelobes  of  high  magnitudes 
are  present  in  the  search  area.  This  difference  may  be  explained  as  follows.  For 
an  isospeed  environment,  no  ray-bending  will  take  place,  versus  the 
upward  refracting  rays  for  a  positive  sound  speed  gradient  environment.  If  the 
rays  do  not  bend,  an  array  can  intercept  multipaths  that  have  traveled  through 
any  depth  of  the  water  column.  For  upward  refracting  rays,  an  array  can  only 
intercept  multipaths  that  have  traveled  through  the  shallower  part  of  the  water 
column.  This  implies  that  for  the  shallow  array  more  multipath  interception  can 
take  place  when  the  rays  do  not  bend,  e.g.,  in  the  isospeed  environment,  and  the 
localization  peak  is  more  unique.  In  the  upward  refracting  case,  the  results  are 
less  unique,  giving  rise  to  more  sidelobe  structure. 

For  the  same  environmental  conditions  for  both  the  transient  and  the 
replicas.  Figs.  4.6  through  4.8  display  the  results  for  the  deep  array  for  all  three 
algorithms.  For  the  IFACM  and  the  KzFACM  algorithm  the  ‘zoomed-in’  ambiguity 
surfaces  are  also  shown. 

Regarding  these  figures  the  following  observations  can  be  made. 
Concerning  the  full  range  search  area,  the  maximum  value  of  the  autocorrelation 
matching  function  for  all  three  algorithms  occurs  at  the  correct  source  location. 
Also,  all  three  results  show  a  very  good  match  and  a  close  resemblance.  The 
IFACM  and  KzFACM  results  show  only  small  sidelobes  of  medium-high 
magnitude  directly  next  to  the  footprint.  The  ZFACM  result  shows  at  this  same 
location  only  a  slightly  larger  sidelobe  with  similar  value.  The  footprints  are 
roughly  equal  in  size.  Comparing  the  ‘zoomed-in’  results  for  IFACM  and 
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Figure  4.6;  For  Case  2,  the  ambiguity  surface  for  the  IFACM  algorithm  for  the  deep  array  with  no 
environmental  mismatch:  (a)  the  full  range  search  area;  (b)  the  ‘zoomed-in’  version. 
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Figure  4.7:  For  Case  2,  the  ambiguity  surface  for  the  KzFACM  algorithm  for  the  deep  array  with 
no  environmental  mismatch:  (a)  the  full  range  search  area;  (b)  the  ‘zoomed-in’  ambiguity  surface. 
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Figure  4.8:  For  Case  2,  the  full  range  ambiguity  surface  for  the  ZFACM  algorithm  for  the  deep 
array  with  no  environmental  mismatch. 


KzFACM  shows  the  same  footprint  size,  although  KzFACM  has  more  areas  with 
meidium  magnitude  sidelobes  surrounding  the  footprint. 

Comparing  the  full  range  search  area  results  for  the  deep  array  with  the 
shallow  array  results  for  Case  2  (Fig.  4.5),  the  following  can  be  noted.  In  both 
cases  the  maximum  occurs  at  the  correct  source  location.  For  all  three 
algorithms,  the  match  for  the  deep  array  is  better  than  for  the  shallow  array,  since 
fewer  sidelobes  are  showing  in  the  search  area.  The  footprint  size  is 
approximately  the  same  for  all  algorithms  for  both  arrays.  This  difference  is 
consistent  with  previous  comments.  In  an  environment  with  a  positive  sound 
speed  gradient,  more  upward  refracting  rays  can  be  intercepted  by  a  deep  array 
than  by  a  shallow  array,  providing  a  more  unique  localization. 
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At  this  point  an  important  observation  can  been  made.  From  the  previous 
thesis  work,  a  major  result  was  that  FACM  was  a  more  robust  algorithm  than 
TACM.  This  means  that  for  these  algorithms  the  actual  envelope  matching  in  the 
time-domain  leads  to  a  better  result.  Using  this  same  argument,  the  algorithms 
developed  for  this  thesis  can  be  compared.  For  all  algorithms  one  parameter  of 
the  envelope  matching  is  the  time-domain.  The  other  parameter  is  either  in  the 
depth-domain  (for  KzFACM),  or  in  the  wavenumber-domain  (for  ZFACM).  For  all 
three  cases  the  resemblance  of  the  results  is  remarkable.  Sidelobes  occur  at 
nearly  the  same  location.  Any  differences  in  results  mainly  concern  sidelobe 
magnitudes  and  relative  peak-to-sidelobe  ratios.  This  suggests  that  for  these 
autocorrelation  matching  algorithms,  the  time-domain  component  of  the 
envelopes  is  the  dominant  factor. 

The  second  situation  for  Case  2  contains  an  environmental  mismatch  in 
the  sound  speed  profile.  The  exact  bottom  profile  and  all  other  input  parameters 
are  known. 

In  a  real  life  scenario,  the  received  transient  signal  would  have  the  correct 
sound  speed  profile  and  the  predictions  (the  replicas)  would  have  the  incorrect 
one,  resulting  in  an  environmental  mismatch.  In  order  to  reduce  the 
computational  load,  the  mismatch  is  implemented  in  the  transient  signal,  leaving 
the  replicas  unchanged.  While  this  is  not  exactly  a  reciprocal  change  for  the 
ambiguity  surface,  the  level  of  performance  degradation  is  expected  to  be  the 
same. 

The  sound  speed  profile  mismatch  is  an  assumed  isospeed  profile  of  1498 
m/s  versus  the  actual  positive  sound  speed  gradient  of  1497  m/s  at  the  surface 
and  1499  m/s  at  a  depth  of  100  m.  Thus,  the  mismatch  is  a  depth  average  of  the 
“true”  profile. 

Figures  4.9  through  4. 1 1  show  the  ambiguity  surfaces  for  the  three 
algorithms  for  the  shallow  array.  For  the  IFACM  and  the  KzACM  algorithms,  the 
‘zoomed-in’  versions  are  shown  as  well. 


42 


0 

10 


i 


30 

40 

1 

50 

60 


•r 


0123456789  10 


4.3  4.35  4.4  4.45  4.5  4.55  4.6  4.65  4.7  4.75  4.8 

Range  (Km) 


Figure  4.9:  For  Case  2,  the  ambiguity  surfaces  for  the  IFACM  algorithm  for  the  shallow  array  in  a 
SSP  mismatch  situation:  (a)  the  full  range  search  area;  (b)  the  ‘zoomed-ln’  ambiguity  surface. 
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Figure  4.1 1 :  For  Case  2,  the  full  range  ambiguity  surface  for  the  ZFACM  algorithm  for  the  shallow 
array  in  a  SSP  mismatch  situation. 

Observing  Figs.  4.9  through  4.11  the  following  remarks  can  be  made. 
Concerning  the  full  range  search  area,  for  all  three  algorithms  the  maximum 
value  of  the  autocorrelation  matching  function  does  not  occur  at  the  source 
location.  The  calculated  source  location  is  approximately  40  m  closer  in  range 
and  1 .5  m  shallower  in  depth  than  the  actual  location.  Besides  that,  the  IFACM 
and  KzFACM  results  show  a  very  good  match  and  a  close  resemblance.  Only 
one  small  sidelobe  is  present  in  the  search  area,  and  the  footprint  size  is  very 
small.  The  ZFACM  result  is  a  moderate  match.  A  number  of  sidelobes  of  high- 
and  medium-high  magnitudes  are  present  on  the  ambiguity  surface,  and  the 
footprint  size  is  twice  as  large  in  range  compared  to  the  IFACM  and  KzFACM 
footprints.  Comparing  the  ‘zoomed-in’  results  of  IFACM  and  KzFACM  displays 
more  clearly  the  source  location  mismatch.  Also,  the  footprint  shapes  show 
some  difference.  The  IFACM  footprint  is  larger  in  range,  but  smaller  in  depth. 
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Comparing  the  shallow  array  results  for  Case  2  with  and  without  the  SSP 
mismatch  (Figs.  4.5  and  4.9-1 1 ),  the  following  can  be  noted.  For  all  three 
algorithms,  the  maximum  values  in  the  ambiguity  surfaces  for  the  SSP  mismatch 
situation  are  significantly  lower  (not  above  0.8)  than  without  mismatch.  However, 
for  the  IFACM  and  KzFACM  algorithms,  fewer  sidelobes  of  medium-  low 
magnitude  appear  in  the  SSP  mismatch  situation  than  without  mismatch. 
Therefore,  remarkably,  the  match  appears  to  be  better  with  a  SSP  mismatch  than 
without.  The  footprint  sizes  are  roughly  the  same  for  both  algorithms  in  both 
situations.  The  ZFACM  algorithm  match  is  less  with  the  SSP  mismatch  than 
without  mismatch.  Although  more  ‘shadows’  appear  without  mismatch,  they  look 
less  scattered  and  have  a  smaller  relative  value  than  with  the  SSP  mismatch. 
Also,  the  ZFACM  footprint  is  smaller  without  mismatch  than  with  the  SSP 
mismatch. 

Figures  4.12  through  4.14  show  the  ambiguity  surfaces  for  the  three 
algorithms  for  the  deep  array.  Again,  for  the  IFACM  and  the  KzACM  algorithms, 
the  ‘zoomed-in’  versions  are  shown  also. 

Regarding  these  figures  the  following  observations  can  be  made. 
Concerning  the  full  range  search  area,  for  all  three  algorithms  the  maximum 
value  of  the  autocorrelation  matching  function  clearly  does  not  occur  at  the 
source  location.  The  maximum  value  of  the  autocorrelation  matching  function  is 
approximately  at  a  range  and  depth  of  8.7  km  and  26  m,  respectively,  where  the 
actual  source  location  is  at  a  range  of  4.6  km  and  at  a  depth  of  27  m.  Since 
numerous  sidelobes  are  present  in  the  search  areas,  none  of  them  being  close  to 
the  actual  source  location,  a  discussion  about  footprint  sizes  is  not  very 
meaningful.  The  ‘zoomed-in’  results  of  IFACM  and  KzFACM  show  a  possible 
match  with  a  close  resemblance  between  these  two  results.  Within  the  selected 
window,  the  calculated  source  location  has  no  range  mismatch  but  there  is  a  7  m 
mismatch  (shallower  than  actual)  in  depth.  However,  this  ‘zoomed-in’  result  is 
very  unlikely  to  be  found,  assuming  that  ‘zooming-in’  normally  is  the  next 
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Figure  4.12:  For  Case  2,  the  ambiguity  surfaces  for  the  IFACM  algorithm  for  the  deep  array  in  a 
SSP  mismatch  situation:  (a)  the  full  range  ambiguity  surface;  (b)  the  ‘zoomed-in’  surface. 
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Figure  4.13:  For  Case  2,  the  ambiguity  surfaces  for  the  KzFACM  algorithm  for  the  deep  array  in  a 
SSP  mismatch  situation:  (a)  the  full  range  search  area;  (b)  the  ‘zoomed-in’  ambiguity  surface. 
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Figure  4.14:  For  Case  2,  the  full  range  ambiguity  surface  for  the  ZFACM  algorithm  for  the  deep 
array  in  a  SSP  mismatch  situation. 


stage  after  ‘full  range  searching’.  Initially  observing  the  full  range  search  area,  a 
‘zoomed-in’  window  would  be  placed  around  displayed  maximum  values.  In  this 
case  that  would  mean  that  the  window  would  be  placed  outside  the  actual  source 
location,  therefore  the  correct  source  location  would  never  be  found.  The  fact 
that  a  possible  match  has  been  found  in  the  ‘zoomed-in’  search  area  is  purely 
based  on  coincidence,  and  on  the  settings  of  the  scenarios.  The  moderate 
footprint  size,  approximately  100  m  in  range  and  5  m  in  depth,  emphasizes  the 
misleading  nature  of  these  results. 

The  most  likely  reason  for  this  complete  mismatch  by  all  three  algorithms 
is  the  fact  that  in  this  upward  refracting  environment  the  deep  array  could  be 
located  in  a  shadow  zone.  These  shadow  zones  are  dependent  on  source  range 
and  depth,  and  any  target  located  in  such  a  zone  is  practicably  ‘invisible’. 
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The  very  close  resemblance  between  the  three  results  is  again 
remarkable.  As  stated  earlier,  the  depth-  or  wavenumber-component  is  of  much 
lesser  importance  than  the  time-component  in  the  envelope  matching  technique 
of  these  autocorrelation  matching  algorithms.  Apparently  this  is  even  the  case 
when  the  match  fails  completely. 

The  last  situation  for  Case  2  contains  an  environmental  mismatch  in  the 
bottom  bathymetry.  This  time,  the  sound  speed  profile  and  all  other  input 
parameters  are  known.  The  bottom  bathymetry  mismatch  is  an  assumed  flat 
bottom  at  a  depth  of  1 00  m  versus  the  actual  bottom. 

Again,  in  order  to  reduce  the  computational  load,  the  mismatch  is 
implemented  in  the  transient  signal,  leaving  the  replicas  unchanged.  Figure  4.15 
shows  the  ambiguity  surfaces  for  the  three  algorithms  for  the  shallow  array,  while 
Fig.  4.16  displays  the  results  for  the  deep  array. 

Observing  Fig.  4.15  the  following  remarks  can  be  made.  For  all  three 
algorithms  the  maximum  value  of  the  autocorrelation  matching  function  occurs  at 
a  range  of  4.72  km  and  a  depth  of  27  m.  The  actual  source  location  is  4.6  km 
and  27  m  deep,  resulting  in  a  range  mismatch  of  120  m  and  a  perfect  match  in 
depth.  All  three  algorithm  results  show  a  close  resemblance  and  a  reasonable 
match  with  the  ZFACM  result  having  sidelobes  of  slightly  higher  magnitudes.  All 
results  show  many  sidelobes.  For  the  IFACM  and  the  KzFACM  results  these 
sidelobes  have  medium  magnitudes  in  front  of  the  source  location,  and  appear 
as  medium-high  magnitude  ‘shadows’  behind  the  source  location.  For  the 
ZFACM  result  these  sidelobes  have  medium-high  magnitudes  in  front  of  the 
source  location  and  the  ‘shadows’  behind  the  source  location  are  ‘red-valued’. 
The  footprints  for  the  IFACM  and  KzFACM  results  are  equally  shaped  and  have 
a  moderate  size.  The  footprint  for  the  ZFACM  result,  however,  is  nearly  twice  as 
large  in  depth. 
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Figure  4.15:  For  Case  2,  the  ambiguity  surfaces  for  the  shallow  array  in  a  bottom  bathymetry 
mismatch  situation:  (a)  the  IFACM  algorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM 
algorithm. 
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Comparing  these  shallow  array  results  with  the  shallow  array  results  with 
no  environmental  mismatch  (Fig.  4.5),  the  following  can  be  noted.  First  of  all,  the 
source  location  is  mismatched  in  range  in  the  bottom  bathymetry  mismatch 
situation  versus  a  perfect  location  match  without  mismatch.  Furthermore,  as  one 
might  expect,  without  mismatch  the  results  give  a  better  match  than  the  results 
with  the  bottom  bathymetry  mismatch.  Note  that  the  ambiguity  surfaces  for  all 
three  algorithms  for  both  situations  show  the  same  pattern.  Even  for  different 
mismatch  situations  there  is  resemblance  between  all  the  results.  However,  with 
the  bottom  bathymetry  mismatch  more  sidelobes  and  ‘shadows’  of  high 
magnitudes  are  present.  Also,  the  footprint  size  for  this  mismatch  situation  is  at 
least  twice  as  large. 

Regarding  the  deep  array  results  in  Fig.  4.16  the  following  observations 
can  be  made.  For  all  three  algorithms  the  maximum  value  of  the  autocorrelation 
matching  function  occurs  at  a  range  of  4.72  km  and  a  depth  of  27  m,  where  the 
actual  source  location  is  4.6  km  and  27  m  deep.  So,  like  the  shallow  array,  this 
results  in  a  range  mismatch  of  120  m  and  a  perfect  match  in  depth.  All  three 
algorithms  show  a  remarkably  good  match.  Note  that  dynamic  range  has  been 
adjusted  downwards  to  show  some  sidelobes  in  the  search  area.  Footprint  sizes 
for  all  three  are  very  small,  although  the  ZFACM  result  shows  more  areas  with 
medium  magnitude  around  the  calculated  source  location. 

Comparing  these  deep  array  results  with  the  shallow  array  results  for  the 
same  mismatch  situation  (Fig.  4.15),  the  following  can  be  noted.  For  both  arrays 
the  source  location  mismatch,  which  is  only  a  range  mismatch,  is  the  same.  The 
deep  array  gives  a  much  better  match  for  all  algorithms.  Practically  no  sidelobes 
are  present  on  the  deep  array  ambiguity  surfaces,  where  the  shallow  array 
ambiguity  surfaces  are  full  of  sidelobes  of  medium-high  magnitudes  and 
‘shadows’  of  high  magnitudes.  Footprint  sizes  are  at  least  four  times  larger  for 
the  shallow  array  results. 
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Comparing  the  deep  array  results  with  the  deep  array  results  without 
environmental  mismatch  (Fig.  4.6),  the  following  can  be  observed.  First  of  all, 
with  the  bottom  bathymetry  mismatch  the  source  location  has  a  slight  range 
mismatch.  Furthermore,  in  both  situations  the  results  of  all  three  algorithms  give 
a  very  good  match.  The  ambiguity  surfaces  show  no  sidelobes,  and  the  footprint 
sizes  are  very  small.  Apparently  the  deep  array  in  a  positive  sound  speed 
environment  is  weakly  sensitive  to  bottom  bathymetry  mismatch  situations,  which 
could  have  been  expected.  Still,  both  shallow  and  deep  array  results  did  produce 
a  range  mismatch  in  the  presence  of  the  bathymetry  mismatch,  indicating  some 
sensitivity. 


C.  CASES:  STRONG  NEGATIVE  SSP  GRADIENT 

For  this  case  the  source  was  located  at  a  range  of  7.2  km  from  the  vertical 
line  array  and  at  a  depth  of  31  m.  The  sound  speed  was  range-independent  with 
a  strong  negative,  bi-linear  gradient  (downward  refracting)  sound  speed  profile  of 
1528  m/s  at  the  surface,  1510  m/s  at  a  depth  of  50  m,  and  1489  m/s  at  a  depth  of 
100  m.  The  density  of  water  was  assumed  to  be  1.0  g/cm^.  The  bottom  was 
chosen  to  be  a  ‘bumpy’  bottom  with  the  following  depths:  102  m  at  0  km;  101  m 
at  2.5  km;  99  m  at  5  km;  102  m  at  7.5  km;  and  99  m  at  10  km.  Its  compressional 
sound  speed  was  1604  m/s,  the  density  1.47  g/cm^,  and  the  attenuation  0.12 
dB/km/Hz. 

Like  the  previous  case,  three  different  environmental  situations  are 
considered  for  Case  3.  For  the  first  situation,  the  environment  is  completely 
known  and  the  transient  and  replica  field  functions  have  exactly  the  same  input 
parameters.  In  other  words,  there  is  no  environmental  mismatch.  Figure  4.17 
shows  the  ambiguity  surfaces  for  the  shallow  array  for  all  three  algorithms. 

Figure  4.18  displays  the  results  for  the  deep  array  for  all  three  algorithms. 
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Figure  4,17:  For  Case  3,  the  ambiguity  surfaces  for  the  shallow  array  with  no  environmental 
mismatch:  (a)  the  IFACM  algorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM  algorithm. 
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Figure  4.18:  For  Case  3,  the  ambiguity  surfaces  for  the  deep  array  with  no  environmental 
mismatch;  (a)  the  IFACM  aigorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM  algorithm. 
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Observing  Fig.  4.17  the  following  remarks  can  be  made.  For  all  three 
algorithms  the  maximum  value  of  the  autocorrelation  matching  function  occurs  at 
the  correct  source  location:  7.2  km  and  31  m  deep.  Ail  three  results  show  a 
strong  resemblance.  The  IFACM  and  KzFACM  results  show  a  reasonably  good 
match  with  some  sidelobes  of  medium-high  magnitudes  present  in  the  search 
area.  The  ZFACM  result  is  a  moderate  match  where  many  sidelobes  of  high 
magnitudes  occur.  The  peak-to-sidelobe  ratio  for  this  result  is  very  small  which 
makes  the  source  localization  difficult.  For  all  three  algorithms,  the  shown 
patterns  are  similar  to  the  expected  ray  paths  for  a  negative  sound  speed  profile. 
The  sidelobes  mainly  occur  in  the  upper  part  of  the  water  column  where  the  rays 
will  cross  each  other  more.  The  IFACM  and  KzFACM  footprints  are  the  same  in 
size  and  shape,  while  the  ZFACM  footprint  is  larger. 

Comparing  these  results  with  the  shallow  array  results  of  Case  1  (Figs. 
4.1(b),  4.2(a)  and  4.3(a)),  the  following  can  be  observed.  Although  in  each  case 
the  transient  and  the  replica  signals  were  100%  the  same,  an  isospeed 
environment  (Case  1 )  provides  a  better  match  for  all  three  algorithms  than  a 
strong  negative  SSP  environment  (Case  3).  The  footprint  sizes  are 
approximately  the  same.  However,  in  Case  3  many  sidelobes  of  high 
magnitudes  are  present  on  the  ambiguity  surface.  This  difference  can  be 
explained  as  follows.  First  of  all,  note  that  the  larger  source  range  for  Case  3 
that  could  contribute  to  the  higher  sidelobes.  Second,  the  downward  refracting 
rays  in  Case  3  will  reduce  the  number  of  multipath  interceptions  by  an  array 
compared  to  an  isospeed  environment  with  no  ray  bending.  Less  multipath 
interception  will  result  in  a  worse  match. 

Regarding  Fig.  4.18  the  following  observations  can  be  made.  Very  high 
values  of  the  autocorrelation  matching  function  for  all  three  algorithms  occur  at 
several  different  locations  in  the  search  area.  Making  an  unrealistic  adjustment 
in  the  dynamic  range,  setting  the  lower  limit  only  0.03  less  than  the  maximum 
value,  results  in  a  good  source  location  match  for  the  IFACM  and  KzFACM 
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algorithms.  However,  for  the  ZFACM  algorithm,  this  adjustment  still  results  in  too 
many  sidelobes  of  high  magnitudes  in  the  search  area.  These  outcomes  are  not 
shown  in  this  report  because  of  the  unrealistic  dynamic  range  adjustment.  Since 
the  numerous  locations  of  sidelobes  result  in  a  failure  to  locate  the  actual  source 
location,  a  discussion  about  the  footprint  sizes  is  not  meaningful. 

Comparing  the  results  for  the  deep  array  with  the  shallow  array  results 
(Fig.  4.17),  the  following  can  be  noted.  Even  without  mismatch,  for  the  deep 
array  all  algorithms  fail  to  result  in  a  reasonable  match,  where  for  the  shallow 
array  all  three  algorithms  provide  some  match.  A  possible  reason  for  this  result 
is  that  the  source  could  be  located  in  a  shadow  zone. 

The  second  situation  for  this  case  contains  an  environmental  mismatch  in 
the  sound  speed  profile.  The  exact  bottom  profile  and  all  other  input  parameters 
are  known. 

The  sound  speed  profile  mismatch  only  exists  in  the  lower  part  of  the 
water  column.  The  sound  speed  profile  of  1528  m/s  at  the  surface  and  1510  m/s 
at  a  depth  of  50  m  is  exactly  the  same  for  both  the  actual  and  the  assumed 
profile.  The  actual  profile  continues  with  a  sound  speed  of  1489  m/s  at  a  depth  of 
100  m,  where  the  assumed  profile  continues  with  the  initial  sound  speed  profile 
gradient  resulting  in  a  sound  speed  of  1492  m/s  at  a  depth  of  100  m.  Thus,  the 
actual  sound  speed  profile  gradient  is  more  negative  (more  downward  refracting) 
over  the  second  part  of  the  water  column  than  the  assumed  sound  speed  profile 
gradient. 

Figure  4.1 9  shows  the  ambiguity  surfaces  of  the  three  algorithms  for  the 
shallow  array.  Observing  this  figure  the  following  remarks  can  be  made.  For  all 
three  algorithms  the  maximum  value  of  the  autocorrelation  matching  function 
occurs  very  close  to  the  actual  source  location.  The  calculated  source  depth  is 
for  all  three  algorithms  approximately  4.5  m  deeper  than  the  actual  depth.  The 
calculated  source  range  for  the  IFACM  algorithm  is  80  m  further  than  the  actual 
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Figure  4.19;  For  Case  3,  the  ambiguity  surfaces  for  the  shallow  array  in  a  SSP  mismatch 
situation:  (a)  the  IFACM  algorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM  aigorithm. 
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range.  The  other  two  algorithms  calculated  the  exact  actual  range.  Note  that  for 
exact  location  determination,  larger  views  are  used  than  shown  in  this  report. 

The  IFACM  and  KzFACM  results  show  a  good  match,  and  the  ZFACM  result  is  a 
moderate  match.  The  ZFACM  ambiguity  surface  clearly  has  higher  sidelobes 
than  the  other  two  surfaces.  The  peak-to-sidelobe  ratio  for  ZFACM  is  less  than 
that  for  the  other  two  algorithms,  which  means  difficult  source  localization.  In  all 
cases  the  footprints  have  a  stretched  shape  in  the  direction  of  the  ray  paths.  The 
footprint  size  for  IFACM  and  KzFACM  is  equal  in  range,  while  KzFACM  is  larger 
in  depth.  ZFACM’s  footprint  is  largest  in  range  and  depth. 

Comparing  the  shallow  array  results  in  a  SSP  mismatch  situation  with  the 
shallow  array  results  for  Case  3  with  no  environmental  mismatch  (Fig.  4.17),  the 
following  can  be  noted.  For  both  situations  the  ambiguity  surfaces  show  patterns 
similar  to  the  ray  path  patterns.  The  patterns  are  essentially  the  same  with 
different  amplitude  levels.  A  mismatch  in  source  depth  appears  in  the  SSP 
mismatch  situation  for  all  algorithms.  The  mismatch  in  source  range  as 
calculated  only  by  the  IFACM  algorithm  is  unexpected  and  can  not  be  explained 
at  this  stage.  All  three  results  show  degradation  in  the  SSP  mismatch  situation, 
where  it  also  must  be  noted  that  this  phenomenon  is  smallest  for  the  KzFACM 
algorithm.  In  other  words,  it  seems  that  KzFACM  in  this  case  is  most  robust  for  a 
SSP  mismatch. 

Figure  4.20  represents  the  ambiguity  surfaces  for  the  three  algorithms  for 
the  deep  array.  Regarding  this  figure  the  following  observations  can  be  made. 
For  all  three  algorithms  the  maximum  value  of  the  autocorrelation  matching 
function  occurs  very  close  to  the  actual  source  location.  The  calculated  source 
depth  is  for  all  three  algorithms  approximately  3  m  shallower  than  the  actual 
depth.  The  calculated  source  range  matches  the  actual  range  for  all  three 
algorithms.  Again,  all  three  results  show  a  close  resemblance,  where  the 
sidelobes  in  the  ZFACM  results  have  higher  values.  The  footprint  sizes  are 
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Figure  4.20:  For  Case  3,  the  ambiguity  surfaces  for  the  deep  array  in  a  SSP  mismatch  situation: 
(a)  the  IFACM  algorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM  algorithm. 
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small,  and  approximately  equal  for  all  three  algorithms.  Remarkable  is  the  fact 
that  the  sidelobes  in  the  search  area  mostly  appear  to  be  in  the  lower  part  of  the 
upper  half  of  the  water  column.  At  that  depth,  the  mismatch  in  sound  speed 
profiles  has  not  yet  been  introduced.  However,  downward  refracting  rays  at  that 
depth  have  experienced  most  of  the  SSP  mismatch  from  the  deeper  part  of  the 
water  column. 

Comparing  the  results  with  the  shallow  array  results  for  the  same  SSP 
mismatch  situation  (Fig.  4.19),  the  following  can  be  noted.  The  different  array 
depths  lead  to  different  mismatches  in  the  source  depth.  The  shallow  array 
algorithm  results  calculate  the  source  depth  4.5  m  deeper  than  the  actual  source 
depth.  The  deep  array  algorithm  results  calculate  the  source  depth  3  m 
shallower  than  the  actual  source  depth.  The  deep  array  results  show  fewer 
sidelobes  in  the  ambiguity  surface,  as  well  as  smaller  footprints  for  all  three 
algorithms.  Based  on  these  observations  it  seems  that  with  a  SSP  mismatch 
situation,  the  algorithms  applied  with  the  deep  array  performs  better  than  with  the 
shallow  array.  For  both  arrays  their  order  of  success  is  KzFACM,  IFACM  and 
ZFACM.  The  better  results  of  the  deep  array  can  be  explained  as  follows.  For  a 
negative  sound  speed  gradient,  meaning  downward  refracting  rays,  the  deep 
array  will  intercept  more  multipaths  than  the  shallow  array.  Although  the  deep 
array  is  closer  to  the  mismatch  region  than  the  shallow  array,  apparently  the 
interception  of  more  multipaths  by  the  deep  array  results  in  a  better  match. 

Comparing  the  deep  array  results  for  the  SSP  mismatch  situation  with  the 
deep  array  results  without  environmental  mismatch  (Fig.  4.18),  the  following  can 
be  observed.  The  SSP  mismatch  situation  results  give  a  better  match  of  the 
actual  source  location  than  without  mismatch.  Although  the  source  location  has 
a  small  depth  mismatch  for  the  SSP  mismatch  situation,  the  maximum  in  the 
ambiguity  surface  appears  in  the  vicinity  of  the  actual  source  location.  Without 
mismatch,  many  sidelobes  of  high  magnitudes  are  present,  and  none  of  them  is 
really  close  to  the  actual  source  location.  The  reason  for  the  considerable 
improvement  in  performance  with  the  SSP  mismatch  is  unclear,  however. 

62 


The  last  situation  for  Case  3  contains  an  environmental  mismatch  in  the 
bottom  bathymetry.  The  sound  speed  profile  and  all  other  input  parameters  are 
known.  The  bottom  bathymetry  mismatch  is  an  assumed  flat  bottom  at  a  depth 
of  100  m  versus  the  actual  bottom.  Figures  4.21  and  4.22  show  the  ambiguity 
surfaces  for  the  three  algorithms  for  the  shallow  array  and  the  deep  array, 
respectively. 

Regarding  Fig.  4.21  the  following  observations  can  be  made.  For  all  three 
algorithms  the  maximum  value  of  the  autocorrelation  matching  function  occurs 
somewhat  near  the  actual  source  location.  The  actual  source  range  of  7.2  km  is 
calculated  80  m  further,  the  actual  source  depth  of  31  m  is  calculated  3  m 
deeper.  Besides  that,  all  three  results  show  a  familiar  pattern  similar  to  the  ray 
path  pattern  for  a  negative  sound  speed  profile.  The  displayed  matches  are  fairly 
good.  The  IFACM  and  KzFACM  perform  best,  showing  only  striping  of  medium 
magnitudes  through  the  search  area.  The  ZFACM  result  is  slightly  worse.  The 
striping  in  the  search  area  is  of  higher  magnitudes.  Footprint  sizes  are  small  for 
all  three  algorithms. 

Comparing  the  shallow  array  results  for  a  bottom  bathymetry  mismatch 
with  the  shallow  array  results  for  Case  3  with  no  environmental  mismatch  (Fig. 
4.17),  the  following  can  be  noted.  For  both  situations  the  ambiguity  surfaces 
show  patterns  similar  to  the  ray  path  patterns.  A  mismatch  in  source  location 
appears  in  the  bottom  bathymetry  mismatch  situation  for  all  algorithms,  as 
observed  in  Case  2  with  bottom  bathymetry  mismatch.  All  three  results  show  a 
better  match  with  bottom  bathymetry  mismatch,  especially  the  ZFACM  result. 

Also  note  that  the  footprint  sizes  are  smaller  in  range  with  the  bottom  bathymetry 
mismatch  situation. 
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Figure  4.21:  For  Case  3,  the  ambiguity  surfaces  for  the  shallow  array  in  a  bottom  bathymetry 
mismatch  situation:  (a)  the  IFACM  aigorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM 
algorithm. 
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Observing  Fig.  4.22  the  following  remarks  can  be  made.  For  the  shown 
dynamic  range  for  all  three  algorithms,  the  maximum  value  of  the  autocorrelation 
matching  function  does  not  occur  at  the  actual  source  location.  The  calculated 
source  depth  is  for  all  three  algorithms  approximately  1.5  m  deeper  than  the 
actual  depth.  The  calculated  source  range  matches  the  actual  range  for  all  three 
algorithms.  Note  that  in  order  to  determine  exactly  the  calculated  source  location 
from  the  information  in  Fig.  4.22,  an  adjusted  dynamic  range  must  be  used.  The 
adjustment  meant  changing  the  lower  limit  to  0.95.  This  is  probably  an 
unrealistic  dynamic  range  adjustment  and  therefore  these  ambiguity  surfaces  are 
not  shown  here.  All  three  algorithms  show  a  moderate  match  with  many 
sidelobes  of  medium-high  magnitudes  present  on  the  ambiguity  surfaces. 
Remarkable  is  that  for  this  situation  the  ZFACM  result  shows  the  fewest 
sidelobes.  The  KzFACM  result  shows  a  better  match  than  the  IFACM  result. 

The  footprint  sizes  are  approximately  the  same  for  all  three  algorithms. 

Comparing  these  deep  array  results  with  the  shallow  array  results  for  the 
same  mismatch  situation  (Fig.  4.21),  the  following  can  be  noted.  Both  arrays 
have  a  source  depth  mismatch.  The  shallow  array  has  a  mismatch  of  3  m 
(deeper),  the  deep  array  results  have  a  mismatch  of  1 .5  m  (deeper).  The  deep 
array  results  have  no  mismatch  in  range.  The  shallow  array  results  have  a  range 
mismatch  of  80  m  (further).  The  shallow  array  results  show  less  sidelobes  on  the 
ambiguity  surface,  as  well  as  smaller  footprints  for  all  three  algorithms.  Based  on 
these  observations  it  seems  that  with  a  bottom  bathymetry  mismatch,  the 
algorithms  applied  with  the  deep  array  perform  better  than  with  the  shallow  array. 
For  both  arrays  their  order  of  success  is  KzFACM,  IFACM  and  ZFACM.  Again,  it 
seems  that  more  multipath  interceptions,  even  with  introduced  mismatch,  give 
better  results  than  those  with  fewer  interceptions. 
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Figure  4.22:  For  Case  3,  the  ambiguity  surfaces  for  the  deep  array  in  a  bottom  bathymetry 
mismatch  situation;  (a)  the  IFACM  algorithm;  (b)  the  KzFACM  algorithm;  (c)  the  ZFACM 
algorithm. 
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Comparing  the  deep  array  results  with  the  bottom  bathymetry  mismatch 
and  the  deep  array  results  without  environmental  mismatch  (Fig.  4.18),  the 
following  can  be  observed.  The  bottom  bathymetry  mismatch  situation  results 
give  a  better  match  of  the  actual  source  location  than  without  mismatch. 

Although  the  source  location  has  a  small  depth  mismatch  for  the  bottom 
bathymetry  mismatch  situation,  the  maximum  in  the  ambiguity  surface  appears  in 
the  vicinity  of  the  actual  source  location.  Without  environmental  mismatch,  many 
sidelobes  with  high  magnitudes  are  present  scattered  over  most  of  the  ambiguity 
surface. 


D.  SUMMARY  OF  FINDINGS 

The  previous  three  sections  have  shown  many  different  results  for  the 
three  algorithms.  The  algorithms  have  been  tested  over  two  different  array 
depths  and  three  different  environmental  cases,  of  which  for  two  cases  mismatch 
situations  have  been  introduced.  In  total  42  different  ambiguity  surfaces  have 
been  presented  and  discussed.  In  order  to  gain  more  insight  in  relative 
performance  of  the  algorithms  the  abbreviated  results  are  tabulated  (Tables  4.1 , 
4.2  and  4.3)  in  a  summarizing  fashion.  Since  the  previous  sections  contained  the 
detail  comments,  this  section  will  summarize  only  overall  remarks  on  findings  that 
appeared  in  all  three  environments. 

Results  from  earlier  thesis  work  had  shown  that  the  FACM  algorithm 
performs  better  than  the  TACM  algorithm.  As  explained  in  the  last  two 
paragraphs  in  Chapter  II,  this  suggests  that  the  matching  of  the  squared 
amplitudes  of  the  signals  in  the  time-domain  give  better  results  than  in  the 
frequency-domain.  In  this  thesis,  it  has  been  shown  that  the  driving  factor  of  a 
successful  autocorrelation  matching  algorithm  is  the  time-domain  component. 
This  conclusion  is  based  upon  the  fact  that  the  developed  algorithms  for  this 
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thesis  in  most  occasions  showed  a  close  resemblance.  The  sidelobes  occurred 
at  the  same  locations  and  the  differences  between  the  algorithms  were 
essentially  sidelobe  magnitudes  and  peak-to-sidelobe  ratios.  Also,  it  has  been 
shown  that,  although  at  a  much  smaller  scale,  the  KzFACM  algorithm  performs 
better  than  the  ZFACM  algorithm  in  most  cases.  This  implies  that  the  matching 
of  the  squared  amplitudes  of  the  signals  in  the  depth-  and  time-domain  give 
better  results  than  in  the  wavenumber-  and  time-domain. 


IFACM 

KzFACM 

ZFACM 

Case  1 

Source  location  matched 

Source  location  matched 

Source  location  matched 

Isospeed  SSP, 

Few  sidelobes  (max  80%) 

Few  sidelobes  (max  80%) 

Some  sidelobes  (max  89%) 

Shallow  array, 

No  mismatch 

Footprint:  80  x  3  m 

Footprint:  80  x  3  m 

Footprint:  120  x  6  m 

Case  1 

Source  location  matched 

Source  location  matched 

Source  location  matched 

Isospeed  SSP, 

No  sidelobes 

No  sidelobes 

One  sidelobe  (max  62%) 

Deep  array, 

No  mismatch 

Footprint:  80  x  3  m 

Footprint:  80  x  3  m 

Footprint:  80  x  4.5  m 

Table  4.1 :  The  summarized  results  for  the  IFACM,  KzFACM  and  ZFACM  algorithms  in  the  Case  1 
environment. 


The  fact  that  a  better  match  was  sometimes  achieved  with  than  without 
environmental  mismatch  was  initially  surprising.  Later  thoughts  on  this 
phenomenon  led  to  the  following  explanation.  When  no  mismatch  has  been 
introduced,  the  source  data  and  the  replica  data  are  ideally  100%  equal. 
Furthermore,  the  propagation  can  have  similar  structure  at  other  locations  in  the 
waveguide.  Matching  these  similar  data  sets  can  lead  to  many  high  value 
correlations  that  will  appear  as  sidelobes  of  high  magnitudes.  In  an  ambiguity 
surface  with  these  high  magnitude  sidelobes,  the  actual  source  location  may 
easily  be  disguised.  When  the  mismatch  has  been  introduced,  the  source  data 
and  the  replica  data  can  be  very  different  for  most  of  the  search  area,  except 
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IFACM 

KzFACM 

ZFACM 

Case  2 

Source  location  matched 

Source  location  matched 

Source  location  matched 

Positive  SSP, 

Shallow  array, 

Few  ‘shadows’  (max  88%) 

Few  ‘shadows’  (max  88%) 

Some  ‘shadows’  (max  93%) 

No  mismatch 

Footprint:  80  x  6  m 

Footprint:  80  x  6  m 

Footprint:  80  x  9  m 

Case  2 

Source  location  matched 

Source  location  matched 

Source  location  matched 

I  Positive  SSP, 

I  Deep  array. 

One  sidelobe  (max  85%) 

One  sidelobe  (max  85%) 

One  sidelobe  (max  85%) 

No  mismatch 

Footprint:  40  x  3  m 

Footprint:  40  x  6  m 

Footprint:  40  x  6  m  j 

Case  2 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch  j 

(40  m  closer  X  1.5  m 

(40  m  closer  x  1 .5  m 

(40  m  closer  X  1.5  m  | 

Positive  SSP, 

Shallow  array, 

shallower) 

shallower) 

shallower)  1 

SSP  mismatch 

One  sidelobe  (max  92%) 

One  sidelobe  (max  90%) 

Some  sidelobes  (max  96%)  j 

Footprint:  80  x  3  m 

Footprint:  80  x  3  m 

Footprint:  160x  3  m  1 

Case  2 

No  source  location  match 

No  source  location  match 

No  source  location  match  1 

Positive  SSP, 

(very  large  range  mismatch) 

(very  large  range  mismatch) 

(very  large  range  mismatch)  I 

Deep  array, 

SSP  mismatch 

Many  sidelobes  (max  98%) 

Many  sidelobes  (max  98%) 

Many  sidelobes  (max  98%) 

Footprint:  n.a. 

Footprint:  n.a. 

Footprint:  n.a. 

Case  2 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch  1 

Positive  SSP, 

(120  m  further) 

(120  m  further) 

(120  m  further)  | 

Shallow  array. 

Many  ‘shadows’ 

Many  ‘shadows’ 

Many  ‘shadows’  1 

Bathymetry  mismatch 

(max  91%) 

Footprint:  360  x  6  m 

(max  91%) 

Footprint:  360  x  6  m 

(max  91%),  some  sidelobes  j 
(max  88%) 

Footprint:  360  x  12  m 

Case  2 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch 

Positive  SSP, 

(120  m  further) 

(120  m  further) 

(120  m  further) 

Deep  array. 

Bathymetry  mismatch 

No  sidelobes 

No  sidelobes 

No  sidelobes 

Footprint:  40  x  3  m 

Foofprint:40  x  3  m 

Footprint:  40  x  3  m 

Table  4.2:  The  summarized  results  for  the  IFACM,  KzFACM  and  ZFACM  algorithms  in  the  Case  2 
environment,  and  its  mismatch  situations. 


near  the  source  location.  A  source  location  mismatch  may  be  anticipated  in  the 

presence  of  bottom  bathymetry  mismatch  (D’Spain,  1999).  However,  due  to  the 

overall  difference  between  the  source  data  and  the  replica  data,  any  appearing 
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sidelobes  may  be  smaller  in  magnitude  than  the  calculated  source  location. 
Therefore,  the  resulting  ambiguity  surface  with  mismatch  can  display  a  better 
result  than  without  mismatch. 

In  this  thesis  no  ambiguity  surfaces  of  the  FACM  algorithm  for  a  single 
hydrophone  are  displayed.  However,  the  input  files  for  the  IFACM  algorithm  are 
essentially  these  single  hydrophone  results.  Therefore,  the  author  claims  to  have 
seen  a  large  number  of  FACM  ambiguity  surfaces  and  can  comment  on  the 
single  hydrophone  results  versus  the  IFACM  results.  It  turns  out  that  the 
hydrophone  depth  has  a  large  impact  on  the  success  of  the  FACM  algorithm. 
Within  the  span  of  a  vertical  line  array,  the  autocorrelation  matching  results 
between,  for  example,  hydrophone  one,  eight  and  sixteen  can  differ  significantly. 
The  IFACM  algorithm  is  basically  an  average  over  the  16  single  hydrophone 
results.  Therefore,  the  statement  that  multiple  hydrophone  algorithms  give  better 
results  is  valid.  However,  the  possibility  of  a  better  result  with  a  single 
hydrophone  located  at  exactly  the  optimum  depth  should  be  mentioned  as  well. 

The  different  algorithms  have  different  runtimes.  The  KzFACM  algorithm 
takes  about  five  hours  and  the  ZFACM  algorithm  takes  about  two  hours  to 
calculate  and  display  their  ambiguity  surface.  The  IFACM  algorithm  takes  less 
than  five  minutes  to  calculate  and  display  its  ambiguity  surface,  once  the  single 
hydrophone  FACM  results  have  been  computed.  These  calculations  take 
approximately  one  hour  per  hydrophone.  If  these  calculations  could  be  done  for 
all  hydrophones  simultaneously,  the  total  runtime  for  the  IFACM  algorithm  would 
be  just  over  one  hour.  In  that  case,  the  IFACM  algorithm  would  be  twice  as  fast 
as  the  ZFACM  algorithm  and  five  times  as  fast  as  the  KzFACM  algorithm.  Since 
for  most  cases  the  IFACM  and  the  KzFACM  results  show  a  close  resemblance, 
the  runtime  difference  can  become  an  important  argument  in  a  more  operational 
environment. 
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IFACM 

KzFACM 

ZFACM 

Case  3 

Source  location  matched 

Source  location  matched 

Source  location  matched 

Negative  SSP, 

Shallow  array, 

Some  striping  (max  85%) 

Some  striping  (max  85%) 

More  striping  (max  92%) 

No  mismatch 

Footprint:  80  x  7.5  m 

Footprint:  80  x  7.5  m 

Footprint:  120x  9  m 

Case  3 

No  source  location  match 

No  source  location  match 

No  source  location  match 

Negative  SSP, 

Deep  array, 

Much  striping  (max  98%) 

Much  striping  (max  98%) 

Much  striping  (max  94%) 

No  mismatch 

Footprint:  n.a. 

Footprint:  n.a. 

Footprint:  n.a. 

Case  3 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch 

Negative  SSP, 

Shallow  array, 

(80  m  further  x  4.5  m 
deeper) 

(4.5  m  deeper) 

(4.5  m  deeper) 

SSP  mismatch 

Mod.  striping  (max  89%) 

Mod.  striping  (max  90%) 

Much  striping  (max  96%) 

Footprint:  80  x  9  m 

Footprint:  80  x  15  m 

Footprint:  120  x  15  m 

Case  3 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch 

Negative  SSP, 

(3  m  shallower) 

(3  m  shallower) 

(3  m  shallower) 

Deep  array, 

SSP  mismatch 

Little  striping  (max  90%) 

Little  striping  (max  90%) 

Some  striping  (max  91%) 

Footprint:  80  x  3  m 

Footprint:  80  x  3  m 

Footprint:  80  x  3  m 

Case  3 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch 

Negative  SSP, 

(80  m  further  x  3  m  deeper) 

(80  m  further  x  3  m  deeper) 

(80  m  further  x  3  m  deeper) 

Shallow  array. 

Bathymetry  mismatch 

Some  striping  (max  87%) 

Little  striping  (max  87%) 

More  striping  (max  90%) 

Footprint:  80  x  6  m 

Footprint:  80  x  6  m 

Footprint:  80  x  6  m 

Case  3 

Source  location  mismatch 

Source  location  mismatch 

Source  location  mismatch 

Negative  SSP, 

(1 .5  m  deeper) 

(1.5  m  deeper) 

(1 .5  m  deeper) 

Deep  array. 

Bathymetry  mismatch 

Much  striping  (max  96%) 

Much  striping  (max  96%) 

Some  striping  (max  94%) 

Footprint:  120  x  7.5  m 

Footprint:  120  x  7.5  m 

Footprint:  120  x  7.5  m 

Table  4.3:  The  summarized  results  for  the  IFACM,  KzFACM  and  ZFACM  algorithms  in  the  Case  3 
environment,  and  its  mismatch  situations. 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

In  this  thesis,  several  algorithms  are  developed  to  localize  a  broadband 
transient-like  point  source  with  a  vertical  line  array.  The  algorithms  are  an 
adaptation  of  a  previously  developed  frequency-domain  autocorrelation  matching 
(FACM)  algorithm,  based  on  matched  field  processing  principles.  The  FACM 
algorithm  was  used  to  localize  a  transient  with  a  single  hydrophone.  Its 
performance  has  been  evaluated  under  environmental  mismatch  situations,  and 
was  determined  to  be  the  most  successful  algorithm  compared  to  other  similar 
algorithms.  For  this  thesis,  the  newly  developed  algorithms  have  incorporated 
the  additional  information  from  the  multiple  elements  in  the  array  and  adapted  the 
FACM  algorithm.  The  algorithms  are  depth-  and  frequency-domain 
autocorrelation  matching  (ZFACM),  wavenumber-  and  frequency-domain 
autocorrelation  matching  (KzFACM),  and  the  incoherent  summation  of  the  single 
element  FACM  results  (IFACM).  These  algorithms  have  been  evaluated  for 
various  different  environmental  situations  and  mismatches. 

For  many  of  the  evaluated  situations  it  turned  out  that  all  three  algorithms 
produce  similar  results.  This  was  the  case  for  good  matches,  but  also  for  less 
successful  results.  Since  the  only  difference  between  the  algorithms  was  in  the 
treatment  of  the  multiple  element  information,  it  is  concluded  that  mainly  the  time- 
component  of  the  matching  envelopes  is  responsible  for  the  general  structure  of 
the  ambiguity  surfaces.  The  additional  wavenumber-  or  depth-component, 
completing  the  matching  envelopes  of  the  new  algorithms,  merely  ‘details’  the 
final  results  in  amplitudes  and  peak-to-sidelobe  ratios.  In  this  ‘detailing’  region 
the  KzFACM  algorithm  generally  produced  a  slightly  better  match  than  the 
ZFACM  algorithm.  The  IFACM  results  were  surprisingly  similar  to  the  KzFACM 
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results.  In  all  cases,  the  results  suggested  that  the  use  of  the  additional 
information  from  multiple  elements  does  improve  the  accuracy  and  robustness  of 
the  localizations. 

Also  noticeable  was  the  difference  in  calculation  times.  The  IFACM 
algorithm  was  twice  as  fast  as  the  ZFACM  algorithm  and  five  times  as  fast  as 
KzFACM  algorithm.  This  significant  difference  can  be  of  operational  importance, 
especially  if  the  results  are  very  similar. 

For  a  slightly  advanced  environment,  it  was  observed  that  in  some  cases 
the  results  appeared  to  be  better  with  than  without  an  introduced  mismatch.  A 
possible  explanation  for  this  phenomenon  is  the  fact  that,  when  no  mismatch  is 
introduced,  the  matching  envelopes  have  almost  the  same  values.  This  results 
in  an  ambiguity  surface  with  many  high-valued  sidelobes,  disguising  the  actual 
source  location.  Introducing  a  mismatch  results  in  matching  envelopes  that  are 
more  distinct,  except  at  the  source  location.  This  produces  an  ambiguity  surface 
with  only  a  small  source  location  mismatch  and  low-valued  sidelobes. 

The  following  observations  show  that  the  placement  of  the  source  and  the 
receiver  within  the  environment,  and  the  environment  itself,  have  a  large  impact 
on  the  performance  of  the  algorithms.  One  observation  is  that  there  were 
significant  differences  in  the  ambiguity  surfaces  computed  using  the  single 
hydrophone  FACM  algorithm  over  the  length  of  the  array.  Although  these  results 
were  not  shown  in  this  report,  it  was  clear  that  the  combined  use  of  array  data 
with  the  algorithms  developed  here  improved  the  peak-to-sidelobe  ratios  for  the 
localization  results.  Another  observation  is  that  for  upward  refracting  rays,  the 
shallow  array  gave  better  results  than  the  deep  array.  For  downward  refracting 
rays  the  deep  array  was  generally  more  successful.  The  number  of  multipath 
interceptions  is  the  driving  factor  for  this  success.  In  mismatch  situations,  more 
stable  multipath  interceptions  give  a  better  result.  However,  if  the  receiving  array 
intercepts  a  large  number  of  multipaths  pertubated  by  the  environmental 
mismatch,  the  results  presumably  degrade.  This  brings  up  the  third  observation 
which  is  that  in  two  occasions  the  deep  array  was  not  successful  at  all.  For  both 
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cases  the  rays  were  either  upward  or  downward  bending,  causing  shadow  zones 
at  certain  locations  in  the  search  area.  Probably  for  these  two  cases,  the  source- 
array  geometry  is  such  that  a  shadow  zone  is  involved,  causing  the  complete 
mismatch. 


B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

Applying  the  developed  algorithms  to  real  data  seems  to  be  the  most 
logical  and  interesting  follow-on.  With  the  possibilities  of  the  MMPE  model  the 
real  life  scenario  could  be  simulated  as  well,  and  results  of  both  real  data  and 
synthetic  data  could  be  compared. 

Another  path  for  future  work  could  be  a  continuation  of  the  simulation  runs 
with  the  algorithms  for  a  more  realistic  environment.  After  this  environment  has 
been  established,  the  sensitivity  of  the  algorithms  for  (other)  environmental 
mismatches  can  be  explored  further. 
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